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Abstract. Using a periodic box calculation it is shown that, owing to 
q ' helicity conservation, a large scale field can only develop on a resistive 

. timescale. This behaviour can be reproduced by a mean-field dynamo 

I with a and m quenchings that are equally strong and 'catastrophic'. 
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After the original papers by Vainshtein & Cattaneo (1992) and Cattaneo & 
Vainshtein (1991) about 'catastrophic' a and n t quenching, respectively, there 
has been a continuous debate about the existence and relevance of a-effect dy- 
namos in astrophysics, i.e. in the high R m regime. Their results were supported 
by calculations that did not however have dynamo action (see also Cattaneo & 
■ Hughes 1996). In the present paper we show, using a periodic box calculation 

that, owing to helicity conservation, a large scale field can only develop on a 



resistive timescale, and that this behaviour can be reproduced by a mean-field 
dynamo with a and rjt quenchings that are equally strong and 'catastrophic' in 
the sense that quenching sets in once the large scale field reaches the fraction 



— 1/2 

R m of the equipartition value. 



2. Helicity constraint and catastrophic quenching 

Large scale magnetic field generation has long been associated with helicity (e.g., 
Moffatt 1978, Krause & Radler 1980). However, the large scale field generated 
by turbulence with kinetic helicity generally possess also magnetic helicity. This 
is a conserved quantity and characterizes the linkage of magnetic flux structures 
with themselves. If the magnetic helicity is zero initially within a certain volume, 
it can only change if there is a loss of magnetic helicity of preferentially one sign 
through the boundaries (Blackman & Field 2000, Kleeorin et al. 2000). Magnetic 
helicity can also change resistively, allowing the large scale magnetic helicity 
to grow while the small scale magnetic helicity is being dissipated, but this 
is a slow process if the magnetic Reynolds number is large. Using only helicity 
conservation and the assumptions that the large scale field, B, possesses helicity, 
one can show that (Brandenburg 2000, hereafter B2000) 

(B 2 ) « hi AB B\ \l - exp (- 2nkl (t - t sat )l . (1) 
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Here, t sat = A -1 ln(Sfi n /Si n i) is the time when the field at small and interme- 
diate scales saturates, B Ui \ and B^ n are initial and final field strengths, A is the 
kinematic growth rate of the dynamo, l? cq is the equipartition field strength, 
and kf is the wavenumber of the forcing. 

Equation (1) is rather general and independent of the actual model of field 

amplification. If the field is not fully helical, then Iab = \(A ■ B)\/(B ) will be 
less than l/k\, where k\ = 2-rr/L is the smallest possible wavenumber. Thus, 
the final field amplitude will be reduced. In any case however the time it takes 
to reach this final amplitude scales with the resistive time, (ryfcf) -1 , where rj is 
the microscopic (not the turbulent!) value. 

Equation (1) suggests that, after t = t sat the large scale field grows at first 

linearly, (B 2 ) = 2rjk 2 (t — t sa t), and then slowly approaches -Bfi n . In Fig. 1 we 

2 

compare the evolution of (B ) for Run 3 of B2000 with the result of Eq. (1). In 
this run we have B ini = 2 x 1CT 9 , B &n = 0.342, and A = 0.067, so t sat = 283. 
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Eq. (3) : 
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Figure 1. Evolution of (B ) for Run 3 (solid line) compared with 
Eq. (1) (dashed line). The dotted line gives the result for an a 2 mean- 
field model (discussed below). The inset shows the evolution on a 
linear scale. Note the excellent agreement between simulation and the 
analytic fits during the saturation phase. n rms = 0.18, k\ = 1, fef = 5, 
r] = 0.002. 
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3. How to model this? 

Exactly the same resistive asymptotic behaviour can be reproduced using an 
a 2 -dynamo with simultaneous a and ^-quenching of the form 

"o Vto ,„s 

a = =2 ' Vt = =o ) l 2 J 

1 + olbB / £?f q 1 + VBB 2 /Bl q 

where as = rjB is assumed. Assuming that the magnetic energy density of the 

— 2 

mean field, B , is approximately uniform (which is well obeyed in the simula- 
tions) we can obtain the solution B = B(t) of the a 2 -dynamo equations in the 
form (B2000) 

B 2 /(l-B 2 /Bl n r^=Bl ie ^, (3) 

where A = |ao|^i — ^to^i is the kinematic growth rate, and ?7to = V + Vto- The 
reason the field saturates in spite of simultaneous quenching of a and r/ t is that 
there remains still the microscopic diffusivity 7], which is not quenched, prevent- 
ing therefore indefinite growth. This determines the saturation field strength in 
terms of the quenching parameters and rjg and yields 

*b = Vb = 4i flaV, (4) 



where B 2 q = ^o(pu 2 ) is characteristic of the magnetic energy at small scales. 
To a good approximation the spectral energy in the small scales is close to that 
in the large scales, so kiB^ n m kfB 2 Thus, as ~ \j{j]k\ki). Since A scales 
with u Tms k{ we have «b ~ u rrns /(r]ki), so the quenching coefficient scales with 
the magnetic Reynolds number based on the box scale, not the forcing scale 
as one might have expected. Taking into account factors of order unity, e.g. 
A « 0.07ii rms A;f (where % = 5) and (2tt)~ 1 , we have as = f]B ~ R m /60, where 
Rm = u Tms L/r]. 

We conclude that Eq. (3) provides an excellent fit to the numerical simula- 
tions (Fig. 1) and one might therefore be tempted to extrapolate to astrophysical 
conditions. However, real astrophysical systems have open boundary conditions 
and we now need to know whether this could alleviate the issue of very long 
timescales for the mean magnetic field. 
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